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1.  Summary 


The  objective  of  this  project  is  to  develop  a  gadolinium-based  GaN  semiconductor 
device  for  neutron  detection  in  a  radiation  environment  with  ability  to  reject  gamma-ray 
background.  Aiming  at  this  objective,  several  goals  have  been  accomplished  during  the 
past  years’  research. 

Workable  GaN  Schottky  diode  radiation  detectors  were  successfully  fabricated  and 
tested.  The  first  year  has,  however,  seen  a  significant  effort  towards  fabricating  a  working 
GaN  radiation  detector.  One  such  detector  design  utilized  semi-insulating  GaN  [1], 
sandwiched  between  ohmic  and  Schottky  contacts.  The  Mg  dopant  added  during  the  GaN 
ammonothennal  growth  resulted  in  deep  level  traps,  and  an  alpha-particle  response  was 
not  observed.  A  second  detector  design,  consisting  of  a  Schottky-ohmic  mesa  structure  on 
a  molecular  beam  epitaxy  (MBE)  grown  super-lattice  Gd:GaN,  also  failed  to  show 
radiation  response  due  to  prohibitively  high  leakage  current  [2]. 

The  breakthrough  was  made  in  the  second  year  when  a  Hydride  Vapour  Phase 
Epitaxy  (HVPE)  grown  GaN  substrate  (Kyma  technology  [3])  was  employed  and  the 
sandwich  structure  Schottky  diode  was  fabricated,  on  which  both  alpha  and  neutron 
response  were  successfully  obtained  [4].  Due  to  wafer’s  high  dislocation  density,  these 
detectors  have  a  relatively  high  leakage  current  which  hindered  their  ability  for  radiation 
detection.  Finally,  a  n-type  HVPE  GaN  substrate  characterized  by  low  carrier 
concentration  and  low  dislocation  density  was  obtained,  detectors  fabricated  on  which 
showed  both  alpha  (241Am  source),  beta  (14C  source),  and  neutron  response.  The  high 
signal-to-noise  ratio  guarantees  their  capability  for  neutron  detection  by  employing  Gd  as 
a  neutron  converter. 
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To  further  understand  the  fundamental  material  electronic  properties  of  these 
detectors,  current-voltage  (IV),  capacitance-voltage  (CV),  and  alpha  spectroscopy 
measurement  were  performed  and  basic  parameters  were  derived,  such  as  the  carrier 
concentration,  Schottky  barrier  height,  depletion  region  depth,  and  charge  collection 
efficiency  (CCE)  etc.  In  addition,  numerical  simulations  were  carried  out  by  using 
Sentaurus  TCAD  software  package[5],  which  helped  us  to  gain  a  deep  insight  into  the 
charge  collection  process  within  the  detector. 

Parallel  to  these  researches  and  findings,  the  proposed  Gd  thin-film  covered  twin- 
detector  design  for  neutron  detection  was  first  validated  by  testing  with  Si  detectors.  The 
key  feature  of  the  design  is  using  a  polyethylene  layer  as  the  internal  conversion  electrons 
(ICE)  separator.  Radiation  tests  at  both  low  and  high  neutron  flux  indicated  that  one 
detector  generates  a  combined  signal  induced  by  both  neutrons  and  gamma  rays,  whereas 
the  other  one  produces  a  signal  induced  only  by  gamma  rays.  Thus  a  net  signal  that  was 
introduced  solely  by  neutrons  can  be  obtained  by  subtracting  the  two  signals.  This 
concept  was  further  validated  by  using  GaN  detector  with  6Li  as  a  neutron  converter.  The 
experiment  by  employing  Gd  converter  and  polyethylene  separator  on  GaN  twin-detector 
device  has  also  been  performed,  and  thus  the  ultimately  objective  of  this  project  is 
reached. 

Three  graduate  students  have  dedicated  to  this  project,  and  their  research  efforts 
culminated  the  production  of  two  MS  theses,  and  one  PhD  thesis.  In  Padhraic  Muligan’s 
thesis[6],  a  multi-detector  digital  neutron  depth  profiling  instrument  was  designed  and 
constructed  in  The  Ohio  State  University  Reactor  Lab  (OSURR)  for  neutron  detections, 
in  Jinghui  Wang’s  thesis[7],  the  Schottky  diode  on  GaN  was  successfully  fabricated;  in 
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Dr.  Praneeth  Kandlakunta’s  thesis[8],  the  twin-detector  scheme  was  proposed  and 
evaluated.  In  addition,  seven  peer-reviewed  journal  papers  have  been  published  [4,  9-14] 
and  two  are  currently  under  preparation  [15,  16],  eight  conference  papers  have  published 
in  proceedings[2,  13,  17-22],  and  a  number  of  presentations  have  been  delivered  in 
various  communities.  The  experiences  working  on  this  project  enhanced  students’ 
abilities  to  undertake  scientific  researches  and  have  shaped  the  directions  of  their  further 
academic  career. 
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2.  Device  Fabrication 


2.1  Device  Fabricated  on  Normal  Growth  GaN 

During  the  past  two  years,  GaN  detectors  have  been  fabricated  on  different  types  of 
semiconductors  (semi-insulating,  n-type)  by  employing  various  structures  (mesa, 
sandwich  [4],  and  super-lattice  [2]).  Experiences  learned  from  evaluating  these  detectors 
are  summarized  as,  1)  semi-insulating  GaN  [1]  does  not  show  radiation  response  due  to 
the  high  density  carrier  trapping  centers  caused  by  Mg  dopant;  2)  Gd  doped  supper-lattice 
GaN  is  hindered  by  its  high  leakage  current  due  to  the  large  amount  of  dislocations 
formed  during  growth  process;  3)  the  one  fabricated  on  n-type  GaN  has  shown  both  alpha 
and  neutron  response[l]. 

Figure  1  shows  the  sandwich  structure  radiation  sensors  constructed  on  freestanding 
n-type  GaN  wafers  procured  from  Kyma  Technologies  Inc.  The  1  cm  x  1  cm  wafers  were 
sliced  from  a  bulk  GaN  boule  grown  via  Hydride  Vapor  Phase  Epitaxy  (HVPE),  and 
subsequently  polished  to  a  nominal  475  pm  thickness.  A  full-backside  ohmic  contact 
(Ti/Al/Ni/Au)  was  deposited  on  the  N-face  of  the  wafer  using  electron  beam  evaporation 
and  a  shadow  mask.  Several  circular  Schottky  diodes  (Ni/Au)  with  concentric  guard  rings 
were  then  deposited  on  the  Ga-face  by  way  of  photolithography.  The  wafer  was  mounted 
in  a  Dual  Inline  Package  (DIP)  using  Ag  paste,  and  the  Schottky  diodes  were  wire 
bonded  to  the  DIP  bonding  leads  with  25-pm  Au  wire.  The  device  was  fabricated  in 
Nanotech  West  Laboratory,  The  Ohio  State  University.  The  employed  nanofabrication 
processes  and  the  related  facilities  are  listed  in  Table  1. 
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Figure  1.  GaN  wafer  with  nine  Schottky  diodes  deposited  on  the  Ga-face,  mounted  in  DIP  (left).  Schematic 
of  single  GaN  Schottky  diode  and  guard  ring  with  device  dimensions,  not  to  scale  (right). 


Table  1.  Nanofabrication  processes  for  GaN  radiation  detector  fabrication 


Nanofabrication  process 

Facility 

Photomask  design 

L-Edit  7.12 

Photolithography 

EV  620  Aligner  (ALN02) 

Asher 

Diener  Pico  Oxygen  Plasma  Asher  (ASH01) 

Etching 

Plasma  Therm  770  SLR  ICP-RIE  system  (ETC04) 

Metal  Deposition 

CHA  Solution  System  E-Gun  Evaporator  (EVP03) 

Rapid  Thermal  Annealing 

AG  Associates  Model  410  Rapid  Thermal  Annealer  (RTA01) 

Wire  Bonding 

Kulicke  and  Soffe  3123  Wire  Bonder  (BND03) 

2.2  Device  Fabricated  on  Special  Growth  GaN 

In  order  to  improve  the  performance  of  our  GaN  detectors,  several  batches  of 
detectors  were  fabricated  on  a  specially  grown  HVPE  n-type  substrate  purchased  from 
Kyma  technology  [3],  which  was  characterized  by  its  low  carrier  concentration  (—10 14 
cm'3)  and  low  dislocation  density  [23].  The  parameters  used  for  the  growth  of  this 
material  are  currently  unknown  by  our  group  for  proprietary  reasons.  These  substrates 
were  specifically  requested  after  a  different  group  published  successful  results  using 
similar  material  (hereafter  referred  to  by  its  production  number  1669)  [23].  To  utilize  as 
much  of  the  limited  special  growth  GaN  as  possible,  1  mm  diameter  diodes  with 
concentric  guard  rings  were  deposited  on  the  wafer,  followed  by  deposition  of  0.5  mm 
diameter  diodes  in  the  open  areas,  see  Figure  2. 
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Figure  2.  (Left)  schematic  of  diode  placement  on  special  growth  GaN  wafer  1669.  (Right)  image  of  GaN 
wafer  with  diodes  deposited. 
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3.  Device  Characterization 


3.1  Current-Voltage 

The  current-voltage  (IV)  characteristic  of  the  fabricated  detectors  on  n-type  GaN  was 
performed  by  using  a  Keithley  2400  SourceMeter,  all  of  which  show  considerable 
improvement  in  leakage  current  compared  to  the  first  year  fabricated  devices.  The  diodes 
fabricated  on  the  special  growth  substrate  showed  a  substantial  decrease  in  current 
density  when  reverse  biased.  Additionally,  the  diodes  could  be  biased  to  much  higher 
voltages  without  approaching  the  breakdown  voltage  of  the  material.  The  special  growth 
material  also  produced  more  successful  diodes  (5  total)  than  the  previously  utilized 
material  (only  1).  As  Figure  3  indicates,  the  I-V  measurements  of  one  of  these  devices 
show  leakage  currents  as  low  as  0.1  nA  constantly  through  a  large  reverse  bias  region. 
This  level  of  leakage  current  is  comparable  with  commercially  available  Si  detector  and 
is  capable  of  producing  high  quality  spectrometry. 


Figure  3.  I-V  measurement  of  one  Schottky  diode. 
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3.2  Capacitance-Voltage 

C-V  measurements  using  a  1  MHz  Boonton  7200  capacitance  meter  on  the  n-type 
GaN  wafer  also  showed  a  slightly  lower  carrier  concentration  in  some  diodes  (9.57  x  10l3 
cm’1)  than  the  detectors  fabricated  on  the  standard  HVPE  material  (1.65  x  1016  cm’3).  One 
of  the  diodes  on  the  special  growth  material  indicated  an  increased  Schottky  barrier 
height,  which  could  be  attributed  to  improvements  in  our  surface  cleaning  methods  prior 
to  metallization. 


Table  2.  Diode  Parameters  Extracted  from  C-V  Profiling 


Special  Growth  Wafer 

Standard  Wafer 

Diode  Number 

1669-N7 

1669-05 

1669-07 

2172-F7 

Doping  (cm'3) 

2.06  x  1016 

2.31  x  1016 

9.57  x  1015 

1.65  x  1016 

Schottky  Barrier  Height  (V) 

1.29 

1.04 

1.08 

1.04 

Using  the  parameters  obtained  from  C-V  profiling  of  the  detectors,  the  depletion 
depth  in  the  detector  as  a  function  of  reverse  bias  was  calculated  assuming  an  abrupt 
junction  approximation.  The  results  indicate  two  of  the  diodes  on  the  special  growth 
wafer  have  depletion  depths  quite  similar  to  that  of  the  standard  wafer.  One  of  the  diodes 
on  the  special  growth  material  had  a  slightly  larger  depletion  depth,  2.28  pm  at  -50  V, 
compared  to  a  maximum  of  1.74  pm  for  the  standard  material.  Although  the  gain  in 
depletion  depth  by  using  the  special  growth  material  was  not  considerable,  the  reduction 
in  leakage  current  was  certainly  beneficial  to  device  performance.  This  became  evident 
when  testing  the  radiation  sensitivity  of  the  detectors  using  alpha  and  beta  particle 
sources. 
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Reverse  Bias  (V) 


Figure  4.  Depletion  depth  of  diodes  on  special  growth  and  standard  GaN  wafers  as  a  function  of  reverse 
bias.  Generated  assuming  an  abrupt  junction  approximation  using  values  derived  from  C-V  measurements 
at  room  temperature. 
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4.  Radiation  Response 

4.1  Alpha  Particle  Spectra  from  Normal  Growth  GaN 

The  experimental  setup  of  the  alpha  spectroscopy  system  is  shown  in  Figure  5,  a 
vacuum  chamber  is  used  to  minimize  the  particles’  energy  loss  between  the  source  and 
the  detector,  and  at  the  same  time  its  functions  as  an  electromagnetic  shield  to  reduce 
noise. 


Figure  5.  Alpha  spectroscopy  system:  1,  Earth  ground;  2,  Shielding  and  vacuum  chamber;  3,  Probe  station; 
4,  Ohmic  contact;  5,  Schottky  Contact;  6,  BNC  cable;  7,  Preamplifier;  8,  CAEN  N6724  Digitizer;  9,  CAEN 
Power  Supply;  10,  Power  supply  to  Preamplifier;  11,  High  voltage  supply  to  detector;  12,  Remote  control 
of  power  supply;  13,  Digitizer  connect  to  computer;  14,  Computer. 

The  initial  radiation  sensitivity  of  the  sensor  fabricated  on  the  normal  growth  GaN 
was  first  tested  using  alpha  particles.  The  GaN  detectors,  mounted  in  a  DIP,  were  placed 
1.5  cm  below  a  0.8  pCi  241  Am  source  in  a  vacuum  chamber.  The  detectors  were  reverse 
biased  by  applying  a  positive  potential  to  the  ohmic  contact  and  grounding  the  diode. 
Signals  generated  in  the  detector  were  collected  by  a  90  mV/MeV  charge  sensitive 
preamplifier,  and  digitized  by  a  CAEN  S.p.A.  100  Msample/s,  14-bit  digitizer.  Firmware 
in  the  digitizer  provided  a  trapezoidal  filter  for  the  digitized  signal,  and  was  subsequently 
histogrammed.  The  collection  efficiency  of  the  device  was  investigated  by  sweeping  the 
bias  from  0  to  18  V.  Results  of  the  alpha  spectrum  are  featured  below  in  Figure  6. 
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Figure  6.  Alpha  spectrum  from  241  Am  source  for  a  range  of  bias  voltages.  Inset  figure  is  simple  schematic 
for  detector  testing,  with  1=1.5  cm. 


As  the  figure  indicates,  increasing  the  detector  bias  resulted  in  a  shift  of  the  alpha 
peak  to  higher  energies.  This  is  to  be  expected  as  the  5.48  MeV  alpha  particle  has  a  range 
of  14.4  pm  as  predicted  by  SRIM,  far  beyond  the  ~  1  pm  Space  Charge  Region  (SCR) 
generated  under  these  voltages.  By  increasing  the  bias  voltage,  and  thus  the  SCR,  more 
electron-hole  pairs  are  generated  in  the  SCR  per  incident  particle.  This  results  in  a  larger 
charge  collected  by  the  detector,  and  a  larger  voltage  pulse  from  the  preamp.  The 
broadened  peaks  for  larger  biases  also  indicate  some  charge  fails  to  be  collected  by  the 
detector,  due  to  the  inclusion  of  more  trapping  centers  in  the  larger  SCR. 

4.2  Neutron  Spectra  from  Normal  Growth  GaN 

To  explore  the  feasibility  of  GaN  as  a  neutron  detection  device,  the  detector  was 
tested  at  the  Ohio  State  University  Research  Reactor  (OSURR)  neutron  beam  facility. 
The  beam  facility  is  composed  of  a  neutron  collimator,  to  filter  fast  neutrons  and  shape 
the  beam,  a  beam  shutter,  a  vacuum  chamber,  and  a  sample  mount  for  neutron  sensitive 
material.  The  450  kW  light  water  reactor  provides  a  8.6  x  106  n/cm2-s  neutron  flux  to  the 
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sample  mount  in  the  facility.  The  detector  was  placed  several  centimeters  below  the 
neutron  beam,  facing  a  0.5  mm  thick,  1.5  cm  x  1.5  cm  sample  of  oxidized  natural  Li  (7.5% 
6Li),  as  shown  in  Figure  7.  The  Li20  sample  was  mounted  behind  a 
polytetrafluoroethylene  (PTFE)  sheet  with  a  1  cm  diameter  aperture  centered  in  the 
neutron  beam. 


Figure  7.  Simplified  schematic  of  GaN  neutron  detection  experiment,  not  to  scale.  1  =  10.4  cm. 


Figure  8  indicates  that  the  sensor  fabricated  on  normal  growth  GaN  could  detect 
charged  particles  induced  by  neutron  absorption  in  6Li.  Triton  and  alpha  particles  formed 
a  peak  in  the  GaN  detector  at  rather  low  energies,  which  is  to  be  expected  considering  the 
small  SCR  in  the  detector.  Although  a  significant,  low-energy  background  component 
was  present,  the  neutron  signal  was  substantially  larger  in  the  detector.  While  using  a 
high-flux  neutron  source  for  testing  may  seem  unconventional,  considering  the  fact  that 
the  spectrum  was  obtained  with  a  thermal  neutron  flux  in  the  order  of  106  n/cm2-s,  the 
experimental  setup  also  yielded  a  small  detector  solid  angle,  counteracting  the  high  flux 
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used  to  test  the  device.  Simulations  performed  using  MCNP  indicate  the  solid  angle  for 
the  GaN  detector  in  this  experiment  was  (1.26±0.12)  x  10"4  sr.  As  such,  with  the  count 


rate  achieved  in  this  test,  it  is  reasonable  to  consider  the  use  of  this  device  for  low-level 


field  applications,  such  as  nuclear  materials  detection. 
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Figure  8.  Neutron  induced  signal  detected  by  GaN  detector  and  neutron  beam  background. 


4.3  Alpha  Particle  Spectra  from  Special  Growth  GaN 

Spectra  from  all  working  diodes  on  the  special  growth  wafer  can  be  seen  below  in 


Figure  9.  All  diodes  were  initially  reverse  biased  to  30  V.  Leakage  current  for  all  diodes 


was  between  5  and  10  nA  at  this  voltage.  As  the  figure  shows,  the  smaller  0.5  mm 


diameter  diodes  N5  and  N7  exhibit  sharp  peaks,  with  total  efficiencies  less  than  the  larger 


diodes  on  the  wafer.  The  centroid  of  N7  appears  in  a  higher  channel  than  that  of  N5, 


possibly  indicating  the  depletion  region  of  N7  is  deeper  than  N5,  resulting  in  more  energy 


collected  in  the  detector.  The  difference  in  peak  area  between  detectors  N5  and  N7  (3.65 
cps  for  N5  vs.  2.29  cps  for  N7)  is  primarily  due  to  the  difference  in  source-detector 
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distance  and  orientation  between  acquisitions.  The  same  reasoning  can  explain  the 
discrepancy  in  curve  area  between  detectors  05,  06,  and  07. 

The  two  peaks  appearing  in  the  05  and  06  spectra  were  deemed  artificial  because 
these  peaks  do  not  appear  to  correlate  with  any  detection  signal,  and  could  be  due  to 
noise  or  signal  reflection  due  to  impedance  mismatch  in  the  detector.  Detector  07  seemed 
somewhat  promising  for  beta  detection,  given  its  single  peak  with  a  relatively  high 
centroid. 


30  40 

ADC  Channel 


Figure  9.  Am-241  spectra  from  operational  diodes  on  wafer. 


4.4  Beta  Particle  Spectra  from  Special  Growth  GaN 

Collection  of  spectra  from  beta  particles  was  also  attempted  on  the  new  diodes  since 
the  proposed  neutron  convertor  is  Gd,  which  emits  fast  electrons  upon  neutron  capture. 
Diode  N7  appeared  to  show  a  good  response  to  a  14C  source  (EMax  =  156  keV,  EAve  = 
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49.47  keV),  as  shown  in  Figure  10.  The  peak  centered  about  channel  2  has  a  much  higher 
count  rate  than  the  two  reference  background  spectra.  While  it  may  appear  strange  that 
the  background  peak  at  low  energy  in  Figure  10  do  not  appear  in  the  spectra  for  24 'Am 
above,  it  should  be  noted  that  the  lower  detection  threshold  in  the  digital  pulse  processing 
software  was  lowered  to  accommodate  the  low  energy  beta  particles.  Beta  spectra  in 
other  diodes  were  not  as  promising.  The  results  indicates  that  not  only  did  the 
background  appear  at  a  greater  count  rate  than  with  the  source  presents,  but  the  spectra 
appears  to  decrease  in  intensity  over  subsequent  acquisitions,  indicating  a  stability 


Figure  10.  Beta  spectra  in  detector  N7  when  reverse  biased  at  30  V. 


4.5  Dual  Detector  Operation 

One  of  the  goals  of  this  work  is  to  create  a  neutron  detector  insensitive  to  gamma 
rays.  One  method  of  creating  such  a  device  is  to  subtract  the  signal  of  a  neutron 
insensitive  detector  from  a  similar  neutron  sensitive  detector,  effectively  reducing  gamma 
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induced  events.  Since  our  fabrication  process  creates  multiple  detectors  on  a  single  wafer, 
we  can  use  two  or  more  devices  simultaneously,  without  dicing  the  wafer.  Subjecting  the 
wafer  to  a  dicing  procedure  could  impart  residual  stress  to  the  wafer,  degrading  the 
device’s  performance.  To  achieve  this  goal,  I-V  measurements  were  repeated  while 
biasing  two  detectors  simultaneously.  Figure  1 1  shows  I-V  measurements  taken  from 
both  diodes  simultaneously,  in  comparison  with  I-V  measurements  obtained  from  each 
diode  one  at  a  time.  Results  indicate  that  one  of  the  diodes  (05)  begins  to  break  down  at  a 
much  lower  voltage  (-30  V  vs.  -43  V)  when  a  second  diode  on  the  wafer  is  also  reverse 
biased.  Similarly,  diode  N7  has  a  slightly  higher  leakage  current  for  reverse  bias  greater 
than  -40  V  if  there  is  a  second  diode  also  biased  on  the  wafer.  Results  from  these 
experiments  helped  determine  the  appropriate  voltage  for  each  diode  when  collecting  an 
alpha  spectrum  from  both  detectors  simultaneously. 


Figure  1 1 .  I-V  curve  of  two  diodes  measured  independently  of  one  another,  and  simultaneously. 
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Spectra  were  first  obtained  using  an  24 'Am  alpha  source  with  both  detectors  biases  to 
-25  V.  Figure  12  shows  the  results  of  this  experiment,  with  the  high  energy  edge  of  both 
peaks  quite  close  to  one  another.  This  is  to  be  expected,  as  the  depletion  depth  predicted 


Figure  12.  Alpha  spectra  obtained  from  diode  N7  and  05  simultaneously  biased  to  -25  V. 


The  experiment  was  repeated  (Figure  13)  with  diode  N7  biased  to  -50  V,  as  the 
leakage  current  in  the  diode  at  this  voltage  was  still  acceptably  low.  As  expected, 
increasing  the  bias  voltage  in  N7  shifted  the  peak  centroid  to  a  higher  channel  and 
broadened  the  spectrum  due  to  the  more  collected  charges  and  increased  statistic  errors. 
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Figure  13.  Alpha  spectra  obtained  from  diode  N7  biased  to  -50  V  and  05  biased  to  -25  V  simultaneously. 

The  results  of  these  experiments  appear  promising  towards  our  goal  of  a  gamma 
discriminating  device.  Future  studies  should  include  determining  the  reason  for  double 
peaks  appearing  in  the  spectra  of  several  diodes,  and  creating  a  process  to  coat  the 
detectors  with  a  neutron  sensitive  material.  Before  the  experiment  on  neutron  converter 
coating,  a  numerical  simulation  is  carried  out  to  the  device  for  understanding  its  charge 
collection  process. 
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5.  Device  Simulation 


5.1  Physical  Models 

In  order  to  further  understand  the  charge  collection  process  in  our  GaN  detectors, 
numerical  simulations  were  perfonned  using  Sentaurus  TCAD  software  package[5].  The 
basic  parameters  obtained  from  the  previous  experiments  are  listed  in  Table  3,  which 
were  incorporated  into  the  simulation  code.  To  minimize  the  calculation  effort,  a  lateral 
size  reduced  2D  structure  was  employed  by  taking  into  account  the  area  factor,  the 
TCAD  codes  are  attached  in  Appendix  1-5. 

The  incident  energetic  particle  is  a  5.48  MeV  alpha  particle  (emitted  by  24 'Am 
source)  and  heavy  ion  model  is  employed.  At  time  of  1><10~15  s,  the  particle  incident  upon 
the  top  Schottky  contact  perpendicularly,  which  is  reverse  biased  at  a  certain  voltage.  The 
particle  then  deposit  its  energy  to  liberate  e-h  pairs  along  its  track  in  GaN  substrate,  the 
linear  transferred  energy  was  calculated  by  SRIM  and  a  simplified  piecewise  step 
function  was  used  to  approximate  the  deposited  energy.  The  carriers  generated  vertically 
along  the  track  are  in  a  Gassian  distribution.  For  carrier  mobility,  constant  mobility  and 
high  field  saturation  models  are  chosen.  For  carrier  recombination,  radiative,  Shockley- 
Read-Hall  (SRH),  and  Auger  recombination  models  are  included.  For  carrier  trapping, 
four  hole  traps  and  four  electron  traps  are  added  with  parameters  taken  from  references 
[24,  25]  shown  in  Table  4  and  Table  5,  respectively. 

Table  3.  Basic  parameters  of  the  GaN-based  Schottky  diode 


Parameter 

Value 

Carrier  concentration  (cm  3) 

1.68xl016  (CV) 

Schottky  barrier  height  (eV) 

0.75  (IV) 

Series  resistance  (Ohm) 

10  (IV) 
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Table  4.  Hole  traps  added  in  the  TCAD  simulation  [24] 


Name 

Energy  from  valance  band  (eV) 

Concentration  (cm  3) 

Capture  cross  section  (cm2) 

H2 

0.55 

2.4E14 

2.7E-12 

H3 

0.65 

2.1E14 

1.7E-14 

H4 

0.85 

3.4E14 

1.8E-13 

H5 

1.2 

2.9E15 

4.7E-14 

Table  5.  Table  5.  Electron  traps  added  in  the  TCAD  simulation  [25] 


Name 

Energy  from  Conduction  band  (eV) 

Concentration  (cm  3) 

Capture  cross  section  (cm2) 

LI 

0.17 

1.0E14 

2.3E-17 

L2 

0.52 

9.0E14 

3.3E-16 

L3 

0.53 

5.8E14 

3.3E-17 

L4 

0.94 

5.8E15 

1.0E-16 

Previous  research  confirmed  that  the  collected  charge  mainly  composed  of  two 
components  from  three  different  regions,  i.e.,  the  drift  and  diffusion  carriers  from 
depletion,  funneling,  and  diffusion  regions.  Figure  14. Structure  of  the  simulated  GaN 
Schottky  diode  and  the  three  regions  for  charge  collection  analysis.  The  purpose  of  this 
simulation  was  to  understand  the  charge  collection  process,  especially  the  CCE 
contributions  of  the  carriers  from  the  three  different  regions. 
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Figure  14.  Structure  of  the  simulated  GaN  Schottky  diode  and  the  three  regions  for  charge  collection 
analysis. 


22 


5.2  Comparison  of  Simulation  and  Experiment  Results 

With  these  physics  models  and  parameters  discussed  before,  Figure  15  compares  the 
simulated  and  measured  capacitance-voltage  (C-V)  and  charge  collection  efficiency 
(CCE)  characteristics  of  the  device.  The  CCE  is  calculated  by  integrating  the  total 
transient  current  on  Schottky  side  in  a  period  of  10  ns.  The  well  agreed  simulation  and 
experiment  results  indicate  that  all  the  major  physics  models  are  correctly  incorporated  in 
our  model. 
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Figure  15.  Simulated  and  measured  results,  left:  capacitance-voltage;  and  right:  charge  collection 
efficiency. 


5.3  Transient  Current  Analysis 

Figure  16  shows  the  total  transient  current  of  the  device  when  reverse  biased  at 
different  voltages.  The  results  indicate  that,  higher  bias  voltage  generate  higher  transient 
current  with  faster  collection  time,  which  can  be  explained  by  the  increased  collected 
charge  (the  width  of  the  depletion  region  is  increased)  and  the  increased  carriers’ 
velocities  (the  electric  field  within  the  semiconductor  region  is  increased),  respectively. 
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Figure  16.  Simulated  total  transient  current  at  different  reverse  bias  voltages  on  the  Schottky  side. 

Figure  17  shows  the  induced  transient  currents  of  different  components  when  the 
diode  is  reversed  biased  at  -20  V,  from  which  it  can  be  concluded  that:  1)  On  Schottky 
side,  the  total  transient  current  is  the  evolution  of  both  holes  and  electrons;  however  on 
Ohmic  side,  the  total  transient  current  is  solely  composed  of  electrons;  2)  On  Schottky 
side,  the  hole  component  goes  through  a  fast  transient;  and  3)  the  total  transient  current  is 
mainly  composed  of  a  fast  ramp  and  a  slow  decay  components.  All  of  these 
characteristics  are  well  explained  by  the  high  built-in  electric  field,  the  thin  depletion 
region,  and  the  lifetimes  and  motilities  of  the  carriers.  For  instance,  even  the  mobility  of 
the  hole  is  smaller  (~  10  cm2/V-s),  they  can  be  collected  in  a  very  short  time  due  to  the 
thin  depletion  region  (1  micrometer).  In  addition,  holes  generated  in  the  bulk  region 
cannot  be  collected  due  to  their  short  lifetime  (~1  ns),  thus  only  the  holes  created  in  the 
depletion  region  contribute  to  the  collected  charge.  That  is  why  the  holes  create  a 
superfast  transient  current  without  a  decay  tail.  On  the  contrary,  electrons  generated  in 
the  depletion  region  and  bulk  region  are  both  contribute  to  the  collected  charge,  and 
electrons  in  the  bulk  region  (low  electric  field)  move  much  slower  than  the  ones 
generated  in  the  depletion  region  (high  built-in  electric  field),  all  of  which  make  the 
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electron  transient  current  characterized  by  a  fast  drift  followed  by  a  slow  diffusion 


component. 
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Figure  17.  Simulated  different  components  of  the  transient  current  created  by  a  5.48  MeV  alpha  particle  at 
a  reverse  bias  of  -20  V. 

When  the  diode  is  reverse  biased  at  -5  V,  Figure  18  and  Figure  19  further  visualized 
the  electron  and  hole  current  density  near  the  Schottky  contact  at  different  transient  times 
induced  by  the  5.48  MeV  alpha  particle,  from  which  the  evolution  of  the  funneling  region 
can  be  clearly  seen. 
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Figure  18.  Simulated  electron  current  densities  at  times  of  le-14,  le-13,  le-12,  5.2e-12,  le-1 1,  and  le-10 
second,  the  units  for  both  x-  and  y-axes  are  micrometer. 


Figure  19.  Simulated  hole  current  densities  at  times  of  le-14,  le-13,  le-12,  5.2e-12,  le-1 1,  and  le-10 
second,  the  units  for  both  x-  and  y-axes  are  micrometer. 
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5.4  Charge  Collection  Efficiency 

Further  simulations  indicate  that  the  width  of  the  funneling  region  is  independent  of 
the  applied  voltage,  i.e.,  about  1  micrometer  for  all  bias  conditions  when  the  hole 
transient  current  reaches  its  maximum.  Thus  we  divided  the  GaN  substrate  into  three 
regions:  depletion  region  (depending  on  the  applied  voltage),  funneling  region  (with  a 
width  of  1  micrometer),  and  diffusion  region  (all  the  remaining  GaN  substrate).  By 
sending  the  alpha  particle  into  different  regions  and  subtracting  the  corresponding 
collected  charges,  the  charge  collection  efficiencies  contributed  by  different  components 
can  be  derived  (see  Figure  20).  The  result  indicates  that  for  Schottky  diode  detectors 
fabricated  on  moderately  doped  semiconductors  (~  1 0 1  6  cm"3  in  our  case),  the  drift  carriers 
in  the  depletion  region  is  still  the  major  component  of  the  collected  charges. 


Figure  20.  Contributions  of  different  transient  current  components  to  the  charge  collection  efficiency. 
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6.  Gamma-ray  Discrimination 

6.1  Gamma-ray  Discrimination  Scheme 

Although  the  small  volume  and  low  Z  number  of  a  thin-film  semiconductor  detector 
(e.g.,  Si)  make  it  relatively  insensitive  to  external  gamma  rays,  those  low-energy  external 
gamma  rays  that  fall  within  the  ICE  energy  range  could  lead  to  false  positive  neutron 
detection.  One  such  application  where  this  could  present  an  issue  is  in  special  nuclear 
material  (SNM)  detection,  which  typically  occurs  in  environments  that  are  “neutron 
signal  starved”,  but  “gamma  signal  abundant”  [26].  To  mitigate  the  likelihood  of  a  false 
detection  due  to  gamma  rays,  a  twin-detector  scheme  is  proposed  using  Gd  as  a  neutron 
converter,  and  two  detectors  to  identify  and  rejecting  external  gamma  rays,  as  shown  in 
Figure  21.  Rejection  is  achieved  by  introducing  a  layer  of  Gd  foil  and  another  layer  of  an 
electron  separator  material  (polyethylene)  into  the  composite  detector  scheme.  Detector  1 
is  placed  in  direct  contact  with  a  Gd  foil,  allowing  the  detection  of  ICEs,  as  well  as 
gamma  rays  produced  by  neutron  absorption  in  Gd.  A  second  detector  is  placed  on  the 
rear  of  the  Gd  layer,  but  uses  a  polyethylene  layer  of  appropriate  thickness  to  stop  all  Gd 
generated  ICEs  and  Auger  electrons  from  reaching  detector  2.  In  this  configuration, 
detector  2  is  sensitive  only  to  gamma  rays.  Hence,  detector  1  generates  a  combined  signal 
induced  by  both  neutrons  and  gamma  rays,  whereas  detector  2  produces  a  signal  induced 
only  by  gamma  rays.  Subtracting  the  two  detector  signals  would  yield  a  net  signal 
induced  solely  by  neutrons. 
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Figure  21.  Gamma  ray  rejection  scheme  with  two  semiconductor  detectors,  a  Gd  layer  and  a  polyethylene 
layer. 


Some  limitations  of  the  proposed  method  were  also  identified  during  a  preliminary 
investigation  and  are  presented  below. 

1)  The  direct  interaction  of  external  gamma  rays  with  detectors  1&2  is  supposed  to 
generate  identical  response  in  the  detectors.  In  addition,  K-X  rays  emitted  from 
gamma  ray  activation  of  Gd  induce  identical  response  in  the  detectors. 
Subtraction  of  the  two  detector  signals  effectively  cancels  out  the  identical 
gamma  ray  component  in  the  resultant  signal.  However,  electrons  emitted  due  to 
photoelectric  absorption  in  Gd  penetrate  detector  1  but  are  blocked  by  the 
polyethylene  separator  from  reaching  detector  2.  Thus,  a  residual  gamma  ray 
component  due  to  photoelectrons  still  exists  in  the  resultant  signal  after 
subtraction. 

2)  Unlike  the  case  in  PSD,  event  by  event  separation  of  neutrons  and  gamma  rays  is 
not  achieved  with  the  proposed  method. 


29 


6.2  Experiment  Setup 

The  proposed  gamma  ray  rejection  scheme  is  heavily  dependent  on  the  detection  of 
ICE  signal  for  registering  neutron  events.  Although  a  large  number  of  prompt  gamma 
rays  in  addition  to  characteristic  x-rays  are  released  in  Gd  neutron  reaction,  the  associated 
component  of  the  neutron  signal  is  effectively  cancelled  out  when  subtracting  the  two 
detector  signals.  In  order  to  validate  our  assumption,  ICE  detections  was  first  done  on  Si 
based  detectors,  and  the  experiment  setup  is  schematic  illustrated  in  Figure  22.  The  whole 
setup  is  incorporated  into  the  OSU  research  reactor,  which  can  provide  a  neutron 
radiation  environment  with  high  gamma-ray  dose. 

In  this  experiment,  a  thin  Gd  foil  (1.25  cm  x  1.25  cm  x  0.0025  cm)  was  mounted  on 
the  sample  holder  inside  the  high-vacuum  chamber.  The  multidetector  setup  provided  an 
opportunity  to  acquire  multiple  spectra  during  one  experiment.  Therefore,  two  of  the 
eight  Si  charged  particle  detectors  were  covered  with  350  pm  thick  polyethylene  caps  to 
shield  them  from  the  ICEs,  while  the  remaining  six  detectors  were  left  unshielded.  A 
background  measurement  was  also  perfonned  at  the  same  instrument  settings  by 
removing  the  Gd  foil  from  the  chamber. 
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Figure  22.  Schematic  illustration  of  the  experimental  setup  and  the  DAQ  electronics  used  in  the  ICE 
measurements. 


6.3  Performance  of  Gd  Covered  Si  Detector  at  Low  Reactor  Power 

For  the  first  experiment,  the  reactor  operated  at  5  kWth  and  delivered  a  thermal 
equivalent  neutron  flux  of  about  9.6  x  104  cm  2-s  1  and  a  gamma  dose  rate  of  about  27 
mRad  hr'1  at  the  sample  location.  A  lower  level  threshold  of  23.3  keV  was  set  on  the 
detector  pulses.  The  measured  results  showed  that  the  energy  spectra  of  the  six  uncovered 
(i i.e .,  unshielded  from  ICEs)  Si  detectors  were  identical,  and  the  energy  spectra  of  the  two 
covered  {i.e.,  shielded  from  ICEs)  Si  detectors  were  also  identical  to  each  other. 

Energy  spectra  representative  of  the  two  sets  detectors  obtained  after  background 
subtraction  are  illustrated  in  Figure  23.  A  highly  intense  and  broad  energy  peak  centered 
at  ~71  keV  is  clearly  seen  from  the  bare  detectors;  in  addition,  two  much  less  intense 
peaks  at  131  keV  and  173  keV  are  also  visible.  The  origin  of  the  71  keV  peak  is  due  to 
the  combined  energy  deposition  of  29-88  keV  ICEs  and  characteristic  x-rays  in  the 
detector.  The  peak  broadening  is  attributed  in  part  to  the  limited  energy  resolution  of  the 
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detector,  but  more  to  the  fact  that  ICEs  lose  energy  when  escaping  the  Gd  foil  before 
reaching  the  Si  detector.  On  the  other  hand,  the  energy  spectrum  from  the  covered 
detectors  indicates  that  the  350  pm  thick  polyethylene  layer  completely  blocked  the  ICEs 
while  still  permitting  gamma/x-rays  to  pass.  An  MCNP5  simulation  confirmed  that  the 
350  pm  thick  polyethylene  cap  used  in  this  study  can  absorb  gamma  rays  only  up  to  ~15 
keV. 


Figure  23.  ICE  energy  spectrum  measured  using  Si  detectors  during  neutron  activation  of  a  thin  Gd  foil. 


The  experimental  results  indicate  that  the  energy  spectra  from  the  bare  detectors  is 
attributed  to  mixed  ICE  and  gamma  ray  response,  whereas  that  from  the  polyethylene 
covered  detectors  is  due  only  to  gamma/x-rays.  Thus,  the  neutron  signal  may  be 
effectively  separated  from  the  gamma  ray  background  by  using  the  proposed  gamma-ray 
rejection  scheme  of  Figure  21.  The  results  also  validate  the  hypothesis  that  the  prompt 
gamma  rays  from  Gd*  are  almost  transparent  to  a  small  volume  semiconductor  detector 
(e.g.,  Si),  while  the  ICEs  constitute  the  principal  neutron  signal  in  such  a  detector. 
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6.4  Performance  of  Gd  Covered  Si  Detector  at  High  Reactor  Power 

The  Gd  ICE  measurement  was  repeated  at  a  higher  neutron  flux  by  increasing  the 
reactor  power  to  250  kWth.  The  thermal  equivalent  neutron  flux  and  the  gamma  dose  rate 
were  observed  as  4.8  x  106  cinqs'1  and  1.35  Rad-hr"1  respectively.  A  lower  level 
threshold  of  20.1  keV  was  set  on  the  detector  pulses.  The  energy  spectra  acquired  from 
the  uncovered  and  covered  detectors  and  the  experimental  geometry  in  inset  are  shown  in 
Figure  24.  As  expected,  due  to  higher  neutron  flux,  the  uncovered  detector  spectrum 
indicated  much  greater  energy  resolution  with  the  three  characteristic  ICE  peaks  clearly 
resolved.  An  interesting  feature  is  the  small  peak  at  ~48  keV  on  the  slope  of  the  main 
peak,  which  is  due  to  the  characteristic  x-rays  of  Gd  activated  by  neutrons  as  well  as 
gamma  rays.  This  peak  is  clearly  identified  in  the  energy  spectrum  of  covered  detector 
due  to  the  separation  of  ICEs. 

The  results  also  indicate  that  the  separation  of  gamma  rays  can  be  effective  even  at  a 
much  higher  gamma  dose  rate.  However,  a  fair  judgment  on  the  efficacy  of  the 
discrimination  technique  may  only  be  possible  with  a  measurement  in  a  gamma  abundant 
and  neutron  deficient  radiation  environment. 
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Figure  24.  ICE  energy  spectrum  measured  using  Si  detectors  during  neutron  activation  of  a  thin  Gd  foil  at  a 
higher  neutron  flux  and  gamma  dose  rate. 


6.5  Performance  of  LiF  Thin  Film  Covered  GaN  Twin-detector 

The  feasibility  of  a  GaN  based,  gamma  discriminating  neutron  detector  was  tested  at 
the  OSU  research  reactor  using  a  6LiF:ZnS  thin  film  as  the  neutron  conversion  material. 
Two  detectors  fabricated  on  the  same  GaN  wafer  were  used  to  simultaneously  detect 
neutron  induced  reactions  in  the  6LiF  converter  and  gamma  interactions.  Both  detectors 
were  first  used  to  observe  alpha  particle  emissions  from  239Pu,  verifying  that  both 
detectors  were  operating  properly,  as  shown  in  Figure  25.  One  of  the  detectors  (07) 
displayed  a  centroid  in  a  higher  channel  than  the  second  detector  (05),  due  to  the  fact  that 
detector  07  has  a  thicker  depletion  region  and  collects  more  electron-hole  pairs  generated 
by  the  incident  alpha  particles. 
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239Pu  Spectra  In  Two  GaN  Detectors 


Figure  25.  Alpha  particle  spectrum  from  239Pu  in  GaN  detectors  07  and  05  to  verify  both  detectors  were 
working  properly. 


The  detectors  were  then  covered  by  two  Teflon  apertures  and  a  thin  film  (0.3  mm) 
6LiF:ZnS  scintillation  screen.  One  of  the  two  apertures  supported  the  scintillation  screen, 
while  the  second  completely  covered  one  of  the  detectors  (05),  leaving  the  second 
detector  with  an  unobstructed  path  to  the  LiF:ZnS  film  (Figure  26). 
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Figure  26.  GaN  detectors  mounted  in  dual  inline  package  with  two  Teflon  apertures.  6LiF:ZnS  converter 
film  is  the  light  green  material  on  the  upper  aperture.  The  lower  Teflon  aperture  leaves  detector  07  exposed 
to  the  film,  and  completely  covers  detector  05. 

The  apparatus  depicted  in  Figure  26  was  then  mounted  in  a  vacuum  chamber  at  the 
OSU  Research  Reactor  neutron  beam  facility.  Electrical  connections  were  made  to  the 
GaN  detectors  and  connected  to  a  charge  sensitive  preamp.  Output  from  the  preamp  was 
connected  to  a  14-bit  digitizer  with  a  trapezoidal  filter,  and  then  to  a  PC.  The  detectors 
were  placed  directly  in  the  3  cm  neutron  beam  with  the  reactor  operating  at  250  kW 
thermal  power  (Figure  27). 


Figure  27.  Detector  setup  in  vacuum  chamber  at  OSU  research  reactor  neutron  beam  facility. 
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Results  of  the  experiment  are  depicted  in  Figure  28.  The  two  charged  particles 
emitted  from  6Li,  3H  at  2727  keV  and  4He  at  2055  keV,  were  clearly  discernable  in  the 
spectrum  collected  by  the  uncovered  detector.  However,  since  the  depletion  depth  of  the 
detector  was  1 .6  pm,  much  less  than  the  range  of  these  particles  in  GaN,  only  a  fraction 
of  the  particles’  initial  energy  was  deposited  in  the  active  region  of  the  detector.  Since  the 
stopping  power  of  the  4He  is  greater  than  the  stopping  power  of  3H,  4He  deposited  more 
energy  in  the  depletion  region  than  3H.  By  performing  a  SRIM  calculation,  it  was 
detennined  that  771  keV  was  deposited  by  the  4He  particle  in  the  depletion  layer,  and  135 
keV  was  deposited  by  the  3H.  These  values  were  then  used  to  perform  an  energy 
calibration  on  the  peaks  in  the  uncovered  detector.  The  covered  detector  did  not  exhibit  a 
response  to  charged  particles,  which  is  expected,  but  did  have  a  significant  low  energy 
signal  below  50  keV,  likely  due  to  gamma  rays  in  the  neutron  beam. 
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Neutron  Induced  Spectrum  from  6Li  Scintillation  Screen 


0  200  400  600  800  1000  1200 
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Figure  28.  Observed  spectrum  from  the  uncovered  detector  (red)  and  covered  detector  (blue)  due  to 
charged  particle  emissions  following  6Li  neutron  capture.  Energy  calibration  was  performed  using  SRIM, 
and  fitting  Gaussian  distributions  to  two  peaks  in  the  uncovered  detector. 


7.  Conclusions 

The  goals  accomplished  and  achievement  made  during  the  past  two-years’  research 
are: 

1)  Sandwiched  structure  Schottky  diode  radiation  sensor  was  successfully  fabricated 
on  n-type  freestanding  GaN,  presenting  both  alpha  and  neutron  response; 

2)  By  employing  a  special  growth  GaN  wafer,  the  electric  properties  of  the 
fabricated  detector  was  improved,  the  level  of  leakage  current  is  comparable  with 
commercially  available  Si  detector  and  is  capable  of  producing  high  quality 
spectrometry.  In  addition  to  alpha  particle,  the  detector  also  showed  beta  particle 
response; 
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3)  Numerical  simulations  using  Sentaurus  TCAD  package  was  performed  to 
understand  the  fundamental  material  and  device  properties  at  picosecond  scale, 
the  charge  collection  process  and  transit  properties  within  the  device  were  studied; 

4)  A  twin-detector  scheme  was  designed  and  the  concept  for  gamma-ray 
discrimination  was  validated  by  using  Si  detectors  at  both  low  and  high  neutron 
flux; 

5)  The  GaN  twin-detector  structure  was  fabricated  and  evaluated  by  using  alpha 
particle  and  6Li,  its  capability  for  neutron  detection  was  validated. 

Two-years’  research  has  brought  the  TRL  of  the  proposed  technology  to  3  or  4.  With 
a  complete  set  of  understanding  for  both  GaN  materials  and  devices  that  are  provided 
with  this  study,  it  is  clear  that  GaN  is  capable  of  serving  as  a  base  material  for  neutron 
detection.  The  challenge  is  obtaining  the  higher  purity  materials  to  form  thicker  depletion 
regions  for  full  spectrum  measurement.  One  patent  application  is  filed.  Three  graduate 
students  have  dedicated  to  this  project,  and  their  research  efforts  culminated  the 
production  of  two  MS  thesis,  one  PhD  thesis,  seven  journal  papers,  and  eight  conference 
proceedings. 
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Appendix  1 :  TCAD  code  for  structure  construction 

;  Sandwich  structure  Schottky  diode  based  on  Kyma  freestanding  n-type  GaN 
;  Write  by  Jinghui  Wang,  04-05-2013 

;  ===================  Outer  Boundary  Creation  ========================= 

Define  dimensions 
(define  wsub  20) 

(define  tsub  430) 

(define  tepi  20) 

(define  tepisub  450) 

(define  sleft  0) 

(define  sright  20) 

GaN  Regions 

(sdegeo:create-rectangle  (position  0.0  0.0  0.0)  (position  wsub  tepisub  0.0)  "GaN"  "EpiSubRegion" ) 


;  ===================  Contact  Definition  Placement  ============== 

(sdegeo:define-contact-set  "top_schottky"  4.0  (color:rgb  1.0  0.0  0.0)  "##") 
(sdegeo:define-contact-set  "bot_ohmic"  4.0  (color:rgb  0.0  1.0  0.0)  "//") 

(sdegeo:insert-vertex  (position  sleft  0.0  0.0)) 

(sdegeo:insert-vertex  (position  sright  0.0  0.0)) 

(sdegeo:define-2d-contact  (find-edge-id  (position  5  0.0  0.0))  "top_schottky") 
(sdegeo:define-2d-contact  (find-edge-id  (position  5  tepisub  0.0))  "bot_ohmic") 


;  ==================  Constant  and  Analytical  Profiles  ========================= 

Definitions 

(sdedr:define-constant-profile  "Sub_Dop_Defn"  "PhosphorusActiveConcentration"  1.68e+16) 
Windows 

(sdedr:define-refeval-window  "Sub_Win"  "Rectangle"  (position  0  0  0)  (position  wsub  tepisub  0.0)) 
Placement 

(sdedr:define-constant-profile-placement  "Sub_Place"  "Sub_Dop_Defn"  "Sub_Win") 


;  ===================  Refinements  ====================== 

Definitions 

(sdedr:define-refinement-size  "Sub_Ref_Defn"  4.0  4.0  0.0  0.8  0.8  0.0  ) 
(sdedr:define-refinement-size  "Epi_Ref_Defn"  1.00.40.0  0.20.1  0.0) 


Windows 

(sdedr:define-refinement-window  "Sub_Win"  "Rectangle"  (position  0.0  0.0  0.0)  (position  wsub  tepisub 

0.0) ) 

(sdedr:define-refinement-window  "Epi  Win"  "Rectangle"  (position  0.0  0.0  0.0)  (position  wsub  tepi  0.0) ) 


Placement 

(sdedr:define-refinement-placement  "Sub_Ref_Place"  "Sub_Ref_Defn"  "Sub_Win" ) 
(sdedr:define-refinement-placement  "Epi_Ref_Place"  "Epi_Ref_Defn"  "Epi_Win" ) 
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;  JC:  Add  mesh  refinement  near  Schottky  electrode  for  proper  construction  of  NLM  in  SDevice 
(sdedr:define-refeval-window  "RefWin. Schottky"  "Rectangle"  (position  sleft  0.0  0)  (position  sright  0.025 
0)) 

(sdedr:define-refinement-size  "RefDef. Schottky"  0.05  0.005  0  0.01  0.001  0  ) 
(sdedr:define-refinement-function  "RefDef. Schottky"  "MaxLenlnt"  "GaN"  "All"  0.001  4) 
(sdedr:define-refinement-placement  "RefPlace. Schottky"  "RefDef. Schottky"  "RefWin. Schottky" ) 

;  JC:  Add  mesh  refinement  to  cover  depletion  region  at  most  negative  bias 

(sdedr:define-refeval-window  "RefWin. Depletion"  "Rectangle"  (position  (-  sleft  0)  0.0  0)  (position  (+ 
sright  0)  2  0)) 

(sdedr:define-refinement-size  "RefDef.Depletion"  0.05  0.05  0  0.05  0.05  0  ) 

(sdedr:define-refinement-placement  "RefPlace. Depletion"  "RefDef.Depletion"  "RefWin. Depletion" ) 

;  WJH:Add  mesh  refinement  for  incident  ion  region 

#(sdedr:define-refeval-window  "RefWin. Iontrack"  "Rectangle"  (position  47  0.0  0)  (position  53  25  0)) 
#(sdedr:define-refinement-size  "RefDef.Iontrack"  0.02  0.02  0  0.01  0.01  0  ) 
#(sdedr:define-refinement-placement  "RefPlace. Iontrack"  "RefDef.Iontrack"  "RefWin.Iontrack" ) 


Resolve  the  junction 

#(sdeaxisaligned: set-parameters  "yCuts"  (list  9.99  10.0  10.01)) 

9 

;  Build  Mesh 

(sde:build-mesh  "snmesh"  "  "  "n@node@_msh") 
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Appendix  2:  TCAD  Code  for  I-V  Simulation 

*****jy  simulation-sdevice  code**** 

File  { 

Grid=  "@tdr@" 

Plot—  "@tdrdat@" 

Current=  "@plot@" 

Output=  "@log@" 

Parameter  "@parameter@" 


Electrode  { 

{  Name="top_schottky"  Schottky  Barrier=  0.75  Voltage=  0.0  } 
{  Name="bot_ohmic"  Voltage=  0.0  Resist=10  } 


Physics  { 

Temperature=  300 
AreaFactor=15700 
Recombination  ( 

Radiative 

SRH( 

DopingDependence 

TempDependence 

#  ElectricField(Lifetime=Hurkx  DensityCorrection=None) 

) 

Auger 

#Avalanche  (vanOverstraetendeMan) 

) 

eBarrierTunneling  "NLM" 

Mobility  ( 

ConstantMobility 

#DopingDependence(Masetti) 

HighF  ieldSaturation 


EffectivelntrinsicDensity  (  OldSlotboom  NoFermi ) 

} 

Physics  (Region=  "EpiSubRegion")  { 

Optics  (OpticalGeneration  (SetConstant  (Value=  lei 6))) 

} 

Plot  { 

eDensity  hDensity 

eMobility  hMobility 

eVelocity  hVelocity 

eCurrent  hCurrent 

eQuasiFermi  hQuasiFermi 

egradQuasiFermi  hgradQuasiFermi 

Doping 

DonorC  one  entration 
AcceptorConcentration 
ElectricField 
Potential 
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SpaceCharge 
ConductionB  andEnergy 
ValenceBandEnergy 
BandGap 

RadiativeRecombination 

SRH 

SRHRecombination 

Auger 

AvalancheGeneration 
eAvalanche  hAvalanche 
#OpticalGeneration 
#OpticalIntensity 


Math  { 

NonLocal  "NLM"  (Electrode="top_schottky"  Length=2e-6  Digits=3  EnergyResolution=0.001) 

Extrapolate 
Digits^  7 

ErrEff(electron)=  le-12 
ErrEff(hole)=  le-12 
RHSmax=  le30 
RHSmin=  le-15 
CdensityMin=  l.e-30 
Notdamped=  20 
Iterations=  15 

eMobilityAveraging=  ElementEdge 
hMobilityAveraging=  ElementEdge 
GeometricDistances 
Parameterlnheritance=  Flatten 

} 

Solve  { 

Poisson 

Coupled  (Iterations=100  LineSearchDamping=le-4)  {  Poisson  Electron  } 

Coupled  (Iterations=100  LineSearchDamping=le-4)  {  Poisson  Electron  Hole  } 

#  IV  from  0  to  -20V 

Quasistationary  (Initialstep=  le-3  Increment=  1.1  Maxstep=  0.2  Minstep=  l.e-15 
Goal  {  Name="top_schottky"  Voltage=-20  }) 

{  Coupled  {  Poisson  Electron  Hole  }  } 

#  IV  from  0  to  5V 
NewCurrentFile=  "IV" 

Quasistationary  (Initialstep=  le-3  Increment=  1.1  Maxstep=  0.025  Minstep=  l.e-15 
Goal  {  Name="top_schottky"  Voltage=  2  }) 

{  Coupled  {  Poisson  Electron  Hole  } 

CurrentPlot  (  Time=  (Range=  (0.00  0.50)  Intervals=  20;  Range=  (0.50  0.75)  Intervals=  20; 
Range=  (0.75  1.00)  Intervals=  10)) 

} 

} 
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Appendix  3:  TCAD  Code  for  C-V  Simulation 

*****CV  simulation-sdevice  code**** 

Device  GaNCV  { 

Electrode  { 

{  Name="top_schottky"  Schottky  Barrier=  0.75  Voltage=  0.0  } 

{  Name="bot_ohmic"  Voltage=  0.0  Resist=10  } 

} 

File  { 

Grid=  "@tdr@" 

Plot=  "@tdrdat@" 

Current=  "@plot@" 

Parameter  "@parameter@" 

} 

Physics  { 

Temperature=  300 
AreaFactor=  39250 
Recombination  ( 

Radiative 

SRH( 

DopingDependence 

TempDependence 

#  ElectricField(Fifetime=Hurkx  DensityCorrection=None) 

) 

Auger 

#Avalanche  (vanOverstraetendeMan) 

) 

eBarrierTunneling  "NEM" 

Mobility  ( 

ConstantMobility 

#DopingDependence(Masetti) 

HighF  ieldSaturation 

) 

EffectivelntrinsicDensity  (  OldSlotboom  NoFermi ) 

} 


Physics  (Region^  "EpiSubRegion"){ 

Traps  ( 

(hNeutral  Conc=2.4el4  Level  EnergyMid=0.55  fromValBand  hXSection=2.7e-12  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=2.7el4  Level  EnergyMid=0.65  fromValBand  hXSection=1.7e-14  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=4.4el4  Level  EnergyMid=0.85  fromValBand  hXSection=1.8e-13  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=4.1el5  Level  EnergyMid=1.20  fromValBand  hXSection=4.7e-14  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(eNeutral  Conc=1.0el4  Level  EnergyMid=0.17  fromCondBand  eXSection=2.3e-17  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=9.0el4  Level  EnergyMid=0.52  fromCondBand  eXSection=3.3e-16  eGfactor=l 


ejfactor=0  eCOnstEmissionRate=0) 
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(eNeutral  Conc=5.8el4  Level  EnergyMid=0.53  fromCondBand  eXSection=3.3e-17  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=5.8el4  Level  EnergyMid=0.94  fromCondBand  eXSection=1.0e-16  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

) 

Optics  (OpticalGeneration  (SetConstant  (Value=  lei 2))) 

} 

}  . . 

Plot  { 

eDensity  hDensity 

eMobility  hMobility 

eVelocity  hVelocity 

eCurrent  hCurrent 

#eQuasiFermi  hQuasiFermi 

#egradQuasiFermi  hgradQuasiFermi 

Doping 

DonorC  one  entration 

AcceptorConcentration 

ElectricField 

Potential 

#SpaceCharge 

#ConductionBandEnergy 

ValenceBandEnergy 

#BandGap 

#RadiativeRecombination 

#SRH 

#SRHRecombination 

#Auger 

#  AvalancheGeneration 

#  eAvalanche  hAvalanche 

#  OpticalGeneration 

#  Opticallntensity 
#eB  arrierT  unneling 

#JC:  View  nonlocal  mesh  in  TDR  file 
NonLocal 

} 

Math  { 

NonLocal  "NLM"  (Electrode="top_schottky"  Length=2e-6  Digits=3  EncrgyResolution=0.001) 

Extrapolate 
Digits=  7 

ErrEff(electron)=  le-12 
ErrEff(hole)=  le-12 
RHSmax=  le30 
RHSmin=  le-15 
CdensityMin=  l.e-30 
Notdamped=  20 
Iterations=  15 

eMobility  A  veraging=  ElementEdge 
hMobility  A  veraging=  ElementEdge 
GeometricDistances 
Parameterlnheritance=  Flatten 


49 


Method=  Blocked 
SubMethod=  Pardiso 


File  { 


Output=  "@log@" 
ACExtract=  "@acplot@" 


}  . 

System  { 

GaNCV  cvcurve  (top_schottky=l  bot_ohmic=2) 

Vsource_pset  vs  (1  0)  {dc=0} 

Vsource_pset  vd  (2  0)  {dc=0} 

} 

Solve  { 

Coupled  (Iterations=100  LineSearchDamping=le-4)  {Poisson} 
Coupled  {  Poisson  Electron  } 

Coupled  {  Poisson  Electron  Hole  } 

#  Schottky  from  0  to  -20V 
Quasistationary  ( 

Initialstep=  le-3 
Increment=  1.1 
Maxstep=  0.1 
Minstep=  l.e-15 

Goal  {  Parameter=  vs.dc  Voltage=  -20  } 

){ 

Coupled  {  Poisson  Electron  Hole  } 

} 

Plot(F  ilePrefix=  "n@node@_vn  15") 

#  Schottky  from  -15  to  0V 
NewCurrentFile=  "AC" 

Quasistationary  ( 

Initialstep=  le-5 
Increment=  1.1 
Maxstep=  0.1 
Minstep=  l.e-15 

Goal  {  Parameter=  vs.dc  Voltage=  0  } 

){ 

ACCoupled  ( 

StartFrequency=  le6 
EndFrequency=  le6 
NumberOfPoints=  1 
Decade 
Node(l  2) 

Exclude(vs  vd) 

) 

{Poisson  Electron  Hole} 

Plot  (Time=  (0;  0.1;  0.2;  0.3;  0.4;  0.5;  0.6;  0.7;  0.8;  0.9)  NoOverwrite) 


*End  of  file 
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Appendix  4:  TCAD  Code  for  Transient  Current  Simulation 

TC AD-TRAN SIENT  CURRENT  FILE: 

*****TranSient  current  generated  by  incident  i0n***** 

File  { 

Grid=  "@tdr@" 

Parameter^  "@parameter@" 

Plot—  "@tdrdat@" 

Current=  "@plot@" 

Output=  "@log@" 


} 

Electrode  { 

{  Name="top_schottky"  Schottky  Barrier=  0.75  Voltage=  (0.0  at  -1,  -18  at  0)  } 
{  Name="bot_ohmic"  Voltage=  0.0  Resist=10} 


Physics  { 

Tempera  ture=  300 
AreaFactor=l 

Recombination  ( 

Radiative 

SRH( 

DopingDependence 

TempDependence 

#ElectricField(Lifetime=Hurkx  DensityCorrection=None) 

) 

Auger 

#Avalanche  (vanOverstraetendeMan) 

) 

eBarrierTunneling  "NLM" 

Mobility  ( 

ConstantMobility 

#DopingDependence(Masetti) 

HighF  ieldSaturation(ElectricField) 


EffectivelntrinsicDensity  (  OldSlotboom  NoFermi ) 

HeavyIon( 

Location=(10,  0,  0) 

Direction=(0,  1,  0) 

Time=le-15 

Length=[l  2  5  9  11  12  13  14  15] 

Wt_hi=[0.285  0.3  0.33  0.39  0.48  0.54  0.6  0.51  0.18] 

LET_f=[0.0095  0.01  0.011  0.013  0.016  0.018  0.02  0.017  0.006] 

Gaussian 

PicoColumb 

) 


Physics  (Region=  "EpiSubRegion")  { 
Traps  ( 
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(hNeutral  Conc=2.4el4  Level  EnergyMid=0.55  fromValBand  hXSection=2.7e-12  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=2.7el4  Level  EnergyMid=0.65  fromValBand  hXSection=1.7e-14  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=4.4el4  Level  EnergyMid=0.85  fromValBand  hXSection=1.8e-13  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=4.1el5  Level  EnergyMid=1.20  fromValBand  hXSection=4.7e-14  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(eNeutral  Conc=1.0el4  Level  EnergyMid=0.17  fromCondBand  eXSection=2.3e-17  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=9.0el4  Level  EnergyMid=0.52  fromCondBand  eXSection=3.3e-16  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=5.8el4  Level  EnergyMid=0.53  fromCondBand  eXSection=3.3e-17  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=5.8el4  Level  EnergyMid=0.94  fromCondBand  eXSection=1.0e-16  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 


) 

Optics  (OpticalGeneration  (SetConstant  (Value=  lei 2))) 


Plot  { 

eDensity  hDensity 

eMobility  hMobility 
eVelocity  hVelocity 
eCurrent  hCurrent 
#eQuasiFermi  hQuasiFermi 
#egradQuasiFermi  hgradQuasiFermi 
Doping 

DonorC  one  entration 

AcceptorConcentration 

ElectricField 

Potential 

SpaceCharge 

#ConductionBandEnergy 

#ValenceBandEnergy 

#BandGap 

#RadiativeRecombination 

#SRH 

#SRHRecombination 

#Auger 

#AvalancheGeneration 
#eAvalanchc  hAvalanche 
#OpticalGeneration 
#OpticalIntensity 

HeavylonCharge 

HeavylonChargeDensity 

HTrappedCharge 

} 

Math  { 

NonLocal  "NLM"  (Electrode=  "top  schottky"  Length=  2e-6  Digits=  3  EnergyResolution=  0.001) 
Extrapolate 
Digits^  7 

ErrEff(electron)=  le-12 
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ErrEff(hole)=  le-12 
RHSmax=  le30 
RHSmin=  le-15 
CdensityMin=  l.e-30 
Notdamped=  20 
Iterations^  15 

eMobilityAveraging=  ElementEdge 
hMobilityAveraging=  ElementEdge 
GeometricDistances 
Parameterlnheritance=  Flatten 

NumberOfThreads=  4 

} 

Solve  { 

Coupled  (Iterations=1000  LineSearchDamping=le-4)  {Poisson} 

Coupled  {Poisson  Electron} 

Coupled  {Poisson  Electron  Hole} 

*ramp  Schottky  from  0V  to  -20V 
NewCurrentPrefix-'RampV" 

Transient 

InitialTime=-l 

FinalTime=0 

InitialStep=le-3 

MaxStep=0.1 

MinStep=le-4 

Increment=l.l 

Decrement=2 

AcceptNewtonParameter(ReferenceStep=  le-12) 

) 

{Coupled  {poisson  electron  hole}} 

*maintain  Schottky  at-  20V,  record  the  transient  current  generated  by  incident  ion 
NewCurrentPrefix="Transient_" 

Transient} 

InitialTime=0.0 

FinalTime=2e-9 

InitialStep=  le-15 

MaxStep=5e-l  1 

MinStep=le-18 

Increment=l.l 

Decrements 

AcceptNewtonParameter(ReferenceStep=  le-12) 

) 

{Coupled  {poisson  electron  hole}} 
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Appendix  5:  TCAD  Code  for  Carrier  Concentration  Simulation 

*****Carrier  concentration  at  different  times  after  an  incident  i0n***** 

File  { 

Grid=  "@tdr@" 

Parameter^  "@parameter@" 

Plot—  "@tdrdat@" 

Current=  "@plot@" 

Output=  "@log@" 


Electrode  { 

{  Name="top_schottky"  Schottky  Barrier=  0.75  Voltage=  (0.0  at  -1,  -160  at  0)  } 
{  Name="bot_ohmic"  Voltage=  0.0  } 


Physics  { 

Temperature=  300 
AreaFactor=39250 
Recombination  ( 

Radiative 

SRH( 

DopingDependence 

TempDependence 

#  ElectricField(Lifetime=Hurkx  DensityCorrection=None) 

) 

Auger 

#Avalanche  (Okuto  Eparallel) 

) 

eBarrierTunneling  "NLM" 

Mobility  ( 

ConstantMobility 

#DopingDependence(Masetti) 

HighF  ieldSaturation 


EffectivelntrinsicDensity  (  OldSlotboom  NoFermi ) 

HeavyIon( 

Location=(10,  0,  0) 

Direction=(0,  1,  0) 

Time=le-15 

Length=[l  2  5  9  11  12  13  14  15] 

Wt_hi=[0.285  0.3  0.33  0.39  0.48  0.54  0.6  0.51  0.18] 

LET_f=[0.0095  0.01  0.011  0.013  0.016  0.018  0.02  0.017  0.006] 

Gaussian 

PicoColumb 

) 


Physics  (Region=  "EpiSubRegion")  { 

Traps  ( 

(hNeutral  Conc=2.4el4  Level  EnergyMid=0.55  fromValBand  hXSection=2.7e-12  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 
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(hNeutral  Conc=2.1el4  Level  EnergyMid=0.65  fromValBand  hXSection=1.7e-14  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=3.4el4  Level  EnergyMid=0.85  fromValBand  hXSection=1.8e-13  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(hNeutral  Conc=2.9el5  Level  EnergyMid=1.20  fromValBand  hXSection=4.7e-14  hGfactor=l 
hJfactor=0  hConstEmissionRate=0) 

(eNeutral  Conc=1.0el4  Level  EnergyMid=0.17  fromCondBand  eXSection=2.3e-17  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=9.0el4  Level  EnergyMid=0.52  fromCondBand  eXSection=3.3e-16  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=5.8el4  Level  EnergyMid=0.53  fromCondBand  eXSection=3.3e-17  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

(eNeutral  Conc=5.8el4  Level  EnergyMid=0.94  fromCondBand  eXSection=1.0e-16  eGfactor=l 
ejfactor=0  eCOnstEmissionRate=0) 

) 

Optics  (OpticalGeneration  (SetConstant  (Value=  lei  2))) 

} 


Plot  { 

eDensity  hDensity 

eMobility  hMobility 
eVelocity  hVelocity 
eCurrent  hCurrent 
eQuasiFermi  hQuasiFermi 
egradQuasiFermi  hgradQuasiFermi 
Doping 

DonorC  one  entration 

AcceptorConcentration 

ElectricField 

Potential 

SpaceCharge 

ConductionB  andEnergy 

ValenceBandEnergy 

BandGap 

RadiativeRecombination 

SRH 

SRHRecombination 

Auger 

#AvalancheGeneration 
#eAvalanchc  hAvalanchc 
#OpticalGeneration 
#OpticalIntensity 

HeavylonCharge 

HeavylonChargeDensity 

} 

Math  { 

NonLocal  "NLM"  ( 

Electrode=  "topschottky" 
Length^  2e-6 
Digits=  3 

EnergyResolution=  0.001 

) 

Extrapolate 
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Digits=  7 

ErrEff(electron)=  le-12 
ErrEff(hole)=  le-12 
RHSmax=  le30 
RHSmin=  le-15 
CdensityMin=  l.e-30 
Notdamped=  20 
Iterations=  15 

eMobilityAveraging=  ElementEdge 

*  uses  edge  mobility  instead  of  element  one  for  electron  mobility 
hMobilityAveraging=  ElementEdge 

*  uses  edge  mobility  instead  of  element  one  for  hole  mobility 
GeometricDistances 

*  when  needed,  compute  distance  to  the  interface  instead  of  closest 

*  point  on  the  interface 
Parameterlnheritance=  Flatten 

*  regions  inherit  parameters  from  materials 


Solve  { 

Poisson 

Coupled  {Poisson  Electron} 

Coupled  {Poisson  Electron  Hole} 

NewCurrentPrefix-'RampV" 

Transient 

InitialTime=-l 

FinalTime=0 

InitialStep=le-3 

MaxStep=0.1 

MinStep=le-5 

Increment=l.l 

Decrement=2 

AcceptNewtonParameter(ReferenceStep=  le-12) 

) 

{Coupled  {poisson  electron  hole}} 

NewCurrentPrefix="Transient_" 

Transient! 

InitialTime=0.0 
FinalTime=5e-10 
InitialStep=  le-15 
MaxStep=5e-l  1 
MinStep=  le-15 
Increment=l.l 
Decrement=2 

AcceptNewtonParameter(ReferenceStep=  le-12) 

) 

{Coupled  {poisson  electron  hole} 

Plot  (  FilePrefix="n30_CarrierDensity_"  Time  =  (  le-14;  le-13;  le-12;  5.2e-12;  le-11;  le-10;  le-9) 
NoOverwrite) 

} 
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